Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases.

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

Systems Analysis of Immune Changes after B-cell Depletion in Autoimmune Multiple Sclerosis.

Multiple sclerosis (MS) is a complex genetically mediated autoimmune disease of the central nervous system where anti-CD20-mediated B cell depletion is remarkably effective in the treatment of early disease. While previous studies investigated the effect of B cell depletion on select immune cell subsets using flow cytometry-based methods, the therapeutic impact on patient immune landscape is unknown. In this study, we explored how a therapy-driven " in vivo perturbation " modulates the diverse immune landscape by measuring transcriptomic granularity with single-cell RNA sequencing (scRNAseq). We demonstrate that B cell depletion leads to cell type-specific changes in the abundance and function of CSF macrophages and peripheral blood monocytes. Specifically, a CSF-specific macrophage population with an anti-inflammatory transcriptomic signature and peripheral CD16 + monocytes increased in frequency post-B cell depletion. In addition, we observed increases in TNFα messenger RNA and protein in monocytes post-B cell depletion, consistent with the finding that anti-TNFα treatment exacerbates autoimmune activity in MS. In parallel, B cell depletion also induced changes in peripheral CD4 + T cell populations, including increases in the frequency of TIGIT + regulatory T cells and marked decreases in the frequency of myelin peptide loaded-tetramer binding CD4 + T cells. Collectively, this study provides an exhaustive transcriptomic map of immunological changes, revealing different mechanisms of action contributing to the high efficacy in B cell depletion treatment of MS.

IL-27 produced during acute malaria infection regulates Plasmodium-specific memory CD4+ T cells.

Malaria infection elicits both protective and pathogenic immune responses, and IL-27 is a critical cytokine that regulate effector responses during infection. Here, we identified a critical window of CD4+ T cell responses that is targeted by IL-27. Neutralization of IL-27 during acute infection with Plasmodium chabaudi expanded specific CD4+ T cells, which were maintained at high levels thereafter. In the chronic phase, Plasmodium-specific CD4+ T cells in IL-27-neutralized mice consisted mainly of CD127+ KLRG1- and CD127- KLRG1+ subpopulations that displayed distinct cytokine production, proliferative capacity, and are maintained in a manner independent of active infection. Single-cell RNA-seq analysis revealed that these CD4+ T cell subsets formed independent clusters that express unique Th1-type genes. These IL-27-neutralized mice exhibited enhanced cellular and humoral immune responses and protection. These findings demonstrate that IL-27, which is produced during the acute phase of malaria infection, inhibits the development of unique Th1 memory precursor CD4+ T cells, suggesting potential implications for the development of vaccines and other strategic interventions.

Increased anti-oxidative action compensates for collagen tissue degeneration in vitiligo dermis.

Vitiligo is a common depigmentation disorder characterized by the selective loss of melanocytes. In our daily clinic experience, we noticed that the skin tightness of hypopigmented lesions would be more evident in comparison to that of uninvolved perilesional skin in vitiligo patients. Therefore, we hypothesized that collagen homeostasis might be maintained in vitiligo lesions, irrespective of the substantial excessive oxidative stress that occurs in association with the disease. We found that the expression levels of collagen-related genes and anti-oxidative enzymes were upregulated in vitiligo-derived fibroblasts. Abundant collagenous fibers were observed in the papillary dermis of vitiligo lesions in comparison to uninvolved perilesional skin by electron microscopy. The production of matrix metalloproteinases that degraded collagen fibers was suppressed. The deposition of acrolein adduct protein, which is a product of oxidative stress, was significantly reduced in vitiligo dermis and fibroblasts. As part of the mechanism, we found upregulation of the NRF2 signaling pathway activity, which is an important defense system against oxidative stress. Taken together, we demonstrated that the anti-oxidative action and collagen production were upregulated and that the collagen degeneration was attenuated in vitiligo dermis. These new findings may provide important clues for the maintenance of antioxidant ability in vitiligo lesions.

Myasthenia gravis-specific aberrant neuromuscular gene expression by medullary thymic epithelial cells in thymoma.

Myasthenia gravis (MG) is a neurological disease caused by autoantibodies against neuromuscular-associated proteins. While MG frequently develops in thymoma patients, the etiologic factors for MG are not well understood. Here, by constructing a comprehensive atlas of thymoma using bulk and single-cell RNA-sequencing, we identify ectopic expression of neuromuscular molecules in MG-type thymoma. These molecules are found within a distinct subpopulation of medullary thymic epithelial cells (mTECs), which we name neuromuscular mTECs (nmTECs). MG-thymoma also exhibits microenvironments dedicated to autoantibody production, including ectopic germinal center formation, T follicular helper cell accumulation, and type 2 conventional dendritic cell migration. Cell-cell interaction analysis also predicts the interaction between nmTECs and T/B cells via CXCL12-CXCR4. The enrichment of nmTECs presenting neuromuscular molecules within MG-thymoma is further confirmed immunohistochemically and by cellular composition estimation from the MG-thymoma transcriptome. Altogether, this study suggests that nmTECs have a significant function in MG pathogenesis via ectopic expression of neuromuscular molecules.

Cell-Free DNA Derived From Neutrophils Triggers Type 1 Interferon Signature in Neuromyelitis Optica Spectrum Disorder.

BACKGROUND AND OBJECTIVES: Recently accumulating evidence suggests the pivotal role of type 1 interferon (IFN-1) signature in the pathogenesis of neuromyelitis optica spectrum disorder (NMOSD). However, the mechanism of the initial trigger that augments IFN-1 pathway in the peripheral immune system of NMOSD has yet to be elucidated. METHODS: Clinical samples were obtained from 32 patients with aquaporin-4 antibody-positive NMOSD and 23 healthy subjects. IFN-1 induction in peripheral blood mononuclear cells (PBMCs) by serum-derived cell-free DNA (cfDNA) was assessed in combination with blockades of DNA sensors in vitro. CfDNA fraction was analyzed for DNA methylation profiles by bisulfite sequencing, elucidating the cellular origin of cfDNA. The induction of neutrophil extracellular trap related cell death (NETosis) was further analyzed in NMOSD and control groups, and the efficacy of pharmacologic intervention of NETosis was assessed. RESULTS: Enhanced IFN-1 induction by cfDNA derived from NMOSD was observed in PBMCs with cofactor of LL37 antimicrobial peptide. DNase treatment, cGAS inhibitor, and Toll-like receptor 9 antagonist efficiently inhibited IFN-1 production. DNA methylation pattern of cfDNA in patients with NMOSD demonstrated that the predominant cellular source of cfDNA was neutrophils. Whole blood transcriptome analysis also revealed neutrophil activation in NMOSD. In addition, enhanced NETosis induction was observed with NMOSD-derived sera, and efficient pharmacologic inhibition of NETosis with dipyridamole was observed. DISCUSSION: Our study highlights the previously unrevealed role of cfDNA predominantly released by neutrophil in the induction of IFN-1 signature in NMOSD and further indicate a novel pharmacologic target in NMOSD.

CCR8-targeted specific depletion of clonally expanded Treg cells in tumor tissues evokes potent tumor immunity with long-lasting memory.

Foxp3-expressing CD25+CD4+ regulatory T cells (Tregs) are abundant in tumor tissues. Here, hypothesizing that tumor Tregs would clonally expand after they are activated by tumor-associated antigens to suppress antitumor immune responses, we performed single-cell analysis on tumor Tregs to characterize them by T cell receptor clonotype and gene-expression profiles. We found that multiclonal Tregs present in tumor tissues predominantly expressed the chemokine receptor CCR8. In mice and humans, CCR8+ Tregs constituted 30 to 80% of tumor Tregs in various cancers and less than 10% of Tregs in other tissues, whereas most tumor-infiltrating conventional T cells (Tconvs) were CCR8- CCR8+ tumor Tregs were highly differentiated and functionally stable. Administration of cell-depleting anti-CCR8 monoclonal antibodies (mAbs) indeed selectively eliminated multiclonal tumor Tregs, leading to cure of established tumors in mice. The treatment resulted in the expansion of CD8+ effector Tconvs, including tumor antigen-specific ones, that were more activated and less exhausted than those induced by PD-1 immune checkpoint blockade. Anti-CCR8 mAb treatment also evoked strong secondary immune responses against the same tumor cell line inoculated several months after tumor eradication, indicating that elimination of tumor-reactive multiclonal Tregs was sufficient to induce memory-type tumor-specific effector Tconvs. Despite induction of such potent tumor immunity, anti-CCR8 mAb treatment elicited minimal autoimmunity in mice, contrasting with systemic Treg depletion, which eradicated tumors but induced severe autoimmune disease. Thus, specific removal of clonally expanding Tregs in tumor tissues for a limited period by cell-depleting anti-CCR8 mAb treatment can generate potent tumor immunity with long-lasting memory and without deleterious autoimmunity.

piRNA/PIWI Protein Complex as a Potential Biomarker in Sporadic Amyotrophic Lateral Sclerosis.

The pathological hallmark of the majority of amyotrophic lateral sclerosis (ALS) cases is the mislocalization and aggregation of TAR DNA-binding protein 43 (TDP-43), an RNA-binding protein. Several studies have attributed disease processes of ALS to abnormal RNA metabolism. However, dysregulated biogenesis of RNA, especially non-coding RNA (ncRNA), is poorly understood. To resolve it, RNA-Seq, biochemical, and immunohistochemical analyses were performed on the pyramidal tract of the medulla oblongata of sporadic ALS (sALS) and control postmortem brain samples. Here, we report perturbation of ncRNA biogenesis in PIWI-interacting RNA (piRNA) in several sALS brain samples associated with TDP-43 pathology. In addition, we confirmed the dysregulation of two PIWI homologs, PIWI-like-mediated gene silencing 1 (PIWIL1) and PIWIL4, which bind to piRNAs to regulate their expression. PIWIL1 was mislocalized and co-localized with TDP-43 in motor neurons of sporadic ALS lumbar cords. Our results imply that dysregulation of piRNA, PIWIL1, and PIWIL4 is linked to pathogenesis of ALS. Based on these results, piRNAs and PIWI proteins are potential diagnostic biomarkers and therapeutic targets of ALS.

High cell surface expression and peptide binding affinity of HLA-DQA1*05:03, a susceptible allele of neuromyelitis optica spectrum disorders (NMOSD).

Neuromyelitis optica spectrum disorder (NMOSD) is a relapsing autoimmune disease characterized by the presence of pathogenic autoantibodies, anti-aquaporin 4 (AQP4) antibodies. Recently, HLA-DQA1*05:03 was shown to be significantly associated with NMOSD in a Japanese patient cohort. However, the specific mechanism by which HLA-DQA1*05:03 is associated with the development of NMOSD has yet to be elucidated. In the current study, we revealed that HLA-DQA1*05:03 exhibited significantly higher cell surface expression levels compared to other various DQA1 alleles, and that its expression strongly depended on the amino acid sequence of the α1 domain, with a preference for leucine at position 75. Moreover, in silico analysis indicated that the HLA-DQ encoded by HLA-DQA1*05:03 preferentially presents immunodominant AQP4 peptides, and that the peptide major histocompatibility complexes (pMHCs) are more energetically stable in the presence of HLA-DQA1*05:03 than other HLA-DQA1 alleles. In silico 3D structural models were also applied to investigate the validity of the energetic stability of pMHCs. Taken together, our findings indicate that HLA-DQA1*05:03 possesses a distinct property to play a pathogenic role in the development of NMOSD.

Localization of KRAS downstream target ARL4C to invasive pseudopods accelerates pancreatic cancer cell invasion.

Pancreatic cancer has a high mortality rate due to metastasis. Whereas KRAS is mutated in most pancreatic cancer patients, controlling KRAS or its downstream effectors has not been succeeded clinically. ARL4C is a small G protein whose expression is induced by the Wnt and EGF-RAS pathways. In the present study, we found that ARL4C is frequently overexpressed in pancreatic cancer patients and showed that its localization to invasive pseudopods is required for cancer cell invasion. IQGAP1 was identified as a novel interacting protein for ARL4C. ARL4C recruited IQGAP1 and its downstream effector, MMP14, to invasive pseudopods. Specific localization of ARL4C, IQGAP1, and MMP14 was the active site of invasion, which induced degradation of the extracellular matrix. Moreover, subcutaneously injected antisense oligonucleotide against ARL4C into tumor-bearing mice suppressed metastasis of pancreatic cancer. These results suggest that ARL4C-IQGAP1-MMP14 signaling is activated at invasive pseudopods of pancreatic cancer cells.

Most cases of pancreatic cancer are detected in the later stages when they are difficult to treat and, as a result, survival is low. Over 90% of pancreatic cancers contain genetic changes that increase the activity of a protein called KRAS. This hyperactive KRAS drives cancer growth and progression. Attempts to treat pancreatic cancer using drugs that reduce the activity of KRAS have so far failed. The KRAS protein can accelerate growth in healthy cells as well as in cancer and it does this by activating various other proteins. Drugs that target some of these other proteins could be more effective at treating pancreatic cancer than the drugs that target KRAS. One of these potential targets is called ARL4C. ARL4C is active during fetal development, but it is often not present in adult tissues. Harada et al. investigated whether the protein is important in pancreatic cancer, and what other roles it has in the body, to better understand if it is a good target for cancer treatment. First, Harada et al. used cells grown in the lab to show that ARL4C contributes to the aggressive spread of human pancreatic cancers. Using mice, Harada et al. also showed that blocking the activity of ARL4C in pancreatic cancers helped to slow their progression. Harada et al.'s results suggest that ARL4C could be a good target for new drugs treating pancreatic cancers. Given that this protein does not seem to have important roles in the cells of adults, targeting it is unlikely to have major side effects. Further investigation of ARL4C in more human-like animal models will help to confirm these results.

Anti-AQP4 autoantibodies promote ATP release from astrocytes and induce mechanical pain in rats.

BACKGROUND: Intractable neuropathic pain is a common symptom of neuromyelitis optica spectrum disorder (NMOSD). However, the underlying mechanism of NMOSD pain remains to be elucidated. In this study, we focused on ATP, which is one of the damage-associated molecular patterns, and also a well-recognized molecule involved in peripheral neuropathic pain. METHODS: We assessed the development of pain symptoms by injecting anti-AQP4 recombinant autoantibodies (rAQP4 IgG) into rat spinal cords. We incubated HEK293 cells expressing AQP4 (HEK-AQP4) and rat astrocytes with rAQP4 IgG and assessed the level of ATP in the supernatant. We performed transcriptome analysis of the spinal cords injected with rAQP4 IgG. Pharmacological inhibition was also applied to investigate the involvement of ATP in the development of neuropathic pain in our rat model. The ATP concentration within the cerebrospinal fluid was examined in patients with NMOSD and other neurological diseases. RESULTS: Development of mechanical allodynia was confirmed in rAQP4 IgG-treated rats. AQP4-Ab-mediated extracellular ATP release from astrocytes was observed in vitro, and pharmacological inhibition of ATP receptor reversed mechanical allodynia in the rAQP4 IgG-treated rats. Furthermore, transcriptome analysis revealed elevation of gene expressions related to several ATP receptors including P2rx4 and IL1B in the spinal cord of rAQP4 IgG-treated rats. In patients, CSF ATP concentration was significantly higher in the acute and remission phase of NMOSD than in multiple sclerosis or other neurological disorders. CONCLUSION: Anti-AQP4 antibody was shown to induce the release of extracellular ATP from astrocytes. The ATP-mediated development of mechanical allodynia was also suggested in rats treated with anti-AQP4 antibody. Our study indicates the pivotal role of ATP in the pain mechanism of NMOSD.

Full-length 16S rRNA gene amplicon analysis of human gut microbiota using MinION™ nanopore sequencing confers species-level resolution.

BACKGROUND: Species-level genetic characterization of complex bacterial communities has important clinical applications in both diagnosis and treatment. Amplicon sequencing of the 16S ribosomal RNA (rRNA) gene has proven to be a powerful strategy for the taxonomic classification of bacteria. This study aims to improve the method for full-length 16S rRNA gene analysis using the nanopore long-read sequencer MinION™. We compared it to the conventional short-read sequencing method in both a mock bacterial community and human fecal samples. RESULTS: We modified our existing protocol for full-length 16S rRNA gene amplicon sequencing by MinION™. A new strategy for library construction with an optimized primer set overcame PCR-associated bias and enabled taxonomic classification across a broad range of bacterial species. We compared the performance of full-length and short-read 16S rRNA gene amplicon sequencing for the characterization of human gut microbiota with a complex bacterial composition. The relative abundance of dominant bacterial genera was highly similar between full-length and short-read sequencing. At the species level, MinION™ long-read sequencing had better resolution for discriminating between members of particular taxa such as Bifidobacterium, allowing an accurate representation of the sample bacterial composition. CONCLUSIONS: Our present microbiome study, comparing the discriminatory power of full-length and short-read sequencing, clearly illustrated the analytical advantage of sequencing the full-length 16S rRNA gene.

VIRTUS: a pipeline for comprehensive virus analysis from conventional RNA-seq data.

SUMMARY: The possibility that RNA transcripts from clinical samples contain plenty of virus RNAs has not been pursued actively so far. We here developed a new tool for analyzing virus-transcribed mRNAs, not virus copy numbers, in the data of bulk and single-cell RNA-sequencing of human cells. Our pipeline, named VIRTUS (VIRal Transcript Usage Sensor), was able to detect 762 viruses including herpesviruses, retroviruses and even SARS-CoV-2 (COVID-19), and quantify their transcripts in the sequence data. This tool thus enabled simultaneously detecting infected cells, the composition of multiple viruses within the cell, and the endogenous host-gene expression profile of the cell. This bioinformatics method would be instrumental in addressing the possible effects of covertly infecting viruses on certain diseases and developing new treatments to target such viruses. AVAILABILITY AND IMPLEMENTATION: : VIRTUS is implemented using Common Workflow Language and Docker under a CC-NC license. VIRTUS is freely available at https://github.com/yyoshiaki/VIRTUS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Proenkephalin+ regulatory T cells expanded by ultraviolet B exposure maintain skin homeostasis with a healing function.

Regulatory T (Treg) cells, expressing CD25 (interleukin-2 receptor α chain) and Foxp3 transcription factor, maintain immunological self-tolerance and suppress various immune responses. Here we report a feature of skin Treg cells expanded by ultraviolet B (UVB) exposure. We found that skin Treg cells possessing a healing function are expanded by UVB exposure with the expression of an endogenous opioid precursor, proenkephalin (PENK). Upon UVB exposure, skin Treg cells were expanded with a unique TCR repertoire. Also, they highly expressed a distinctive set of genes enriched in "wound healing involved in inflammatory responses" and the "neuropeptide signaling pathway," as indicated by the high expression of Penk. We found that not only was PENK expression at the protein level detected in the UVB-expanded skin Treg (UVB-skin Treg) cells, but that a PENK-derived neuropeptide, methionine enkephalin (Met-ENK), from Treg cells promoted the outgrowth of epidermal keratinocytes in an ex vivo skin explant assay. Notably, UVB-skin Treg cells also promoted wound healing in an in vivo wound closure assay. In addition, UVB-skin Treg cells produced amphiregulin (AREG), which plays a key role in Treg-mediated tissue repair. Identification of a unique function of PENK+ UVB-skin Treg cells provides a mechanism for maintaining skin homeostasis.

Regulatory T Cell-Specific Epigenomic Region Variants Are a Key Determinant of Susceptibility to Common Autoimmune Diseases.

The contribution of FOXP3-expressing naturally occurring regulatory T (Treg) cells to common polygenic autoimmune diseases remains ambiguous. Here, we characterized genome-wide epigenetic profiles (CpG methylation and histone modifications) of human Treg and conventional T (Tconv) cells in naive and activated states. We found that single-nucleotide polymorphisms (SNPs) associated with common autoimmune diseases were predominantly enriched in CpG demethylated regions (DRs) specifically present in naive Treg cells but much less enriched in activation-induced DRs common in Tconv and Treg cells. Naive Treg cell-specific DRs were largely included in Treg cell-specific super-enhancers and closely associated with transcription and other epigenetic changes in naive and effector Treg cells. Thus, naive Treg cell-specific CpG hypomethylation had a key role in controlling Treg cell-specific gene transcription and epigenetic modification. The results suggest possible contribution of altered function or development of natural Treg cells to the susceptibility to common autoimmune diseases.

Genome-Wide Natural Selection Signatures Are Linked to Genetic Risk of Modern Phenotypes in the Japanese Population.

Elucidation of natural selection signatures and relationships with phenotype spectra is important to understand adaptive evolution of modern humans. Here, we conducted a genome-wide scan of selection signatures of the Japanese population by estimating locus-specific time to the most recent common ancestor using the ascertained sequentially Markovian coalescent (ASMC), from the biobank-based large-scale genome-wide association study data of 170,882 subjects. We identified 29 genetic loci with selection signatures satisfying the genome-wide significance. The signatures were most evident at the alcohol dehydrogenase (ADH) gene cluster locus at 4q23 (PASMC = 2.2 × 10-36), followed by relatively strong selection at the FAM96A (15q22), MYOF (10q23), 13q21, GRIA2 (4q32), and ASAP2 (2p25) loci (PASMC < 1.0 × 10-10). The additional analysis interrogating extended haplotypes (integrated haplotype score) showed robust concordance of the detected signatures, contributing to fine-mapping of the genes, and provided allelic directional insights into selection pressure (e.g., positive selection for ADH1B-Arg48His and HLA-DPB1*04:01). The phenome-wide selection enrichment analysis with the trait-associated variants identified a variety of the modern human phenotypes involved in the adaptation of Japanese. We observed population-specific evidence of enrichment with the alcohol-related phenotypes, anthropometric and biochemical clinical measurements, and immune-related diseases, differently from the findings in Europeans using the UK Biobank resource. Our study demonstrated population-specific features of the selection signatures in Japanese, highlighting a value of the natural selection study using the nation-wide biobank-scale genome and phenotype data.

Innate Myeloid Cell Subset-Specific Gene Expression Patterns in the Human Colon are Altered in Crohn's Disease Patients.

BACKGROUND/AIMS: There is a heterogeneous subset innate myeloid cells, such as macrophages and dendritic cells, in the human intestinal lamina propria. Several studies have demonstrated that these cells contribute to the maintenance of gut homeostasis through the induction of inflammatory responses and tolerance via cell type-specific mechanisms; whereas, disrupted innate immune responses are implicated in the pathogenesis of Crohn's disease (CD). However, the detailed mechanisms by which each innate myeloid subset regulates gut homeostasis and inflammation largely remain unknown. We aimed to clarify the comprehensive gene expression profiles of innate myeloid cell -subsets in the lamina propria from normal human colons (NC) and the inflamed colon sites from patients with Crohn's disease (CDi). METHODS: We performed RNA-sequencing analysis and precise bioinformatics analysis on 3 innate myeloid cell subsets, CD14-CD11c-, CD14-CD11c+, and CD14+CD11c+CD163low cells from NC and CDi. RESULTS: Transcriptional analysis of the 3 subsets from the NC showed distinct gene expression patterns and gene ontology (GO) enrichment analysis revealed the associated innate myeloid subset-specific biological process (BP) terms. In addition, changes in gene expression patterns were observed in innate myeloid subsets from CDi. Furthermore, the core GO-BP terms for the genes upregulated in the innate myeloid cells from CDi were distinct from those found in NC. CONCLUSION: Our data identified the innate myeloid cell subset-specific transcriptomes and the associated enriched GO-BP terms in the NC and found these patterns were altered in CDi.